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Cosmic-ray electrons and positrons propagating in the Galaxy produce diffuse gamma-rays via the 
inverse Compton (IC) process. Tlie low energy target plioton populations with which the cosmic-rays 
interact during propagation are produced by stars, this stellar light being reprocessed by Galactic 
dust. Detailed modelling of the Galactic stellar distribution, dust distribution, and treatment of the 
absorption and scattering of light is therefore required to obtain accurate models for the low energy 
Galactic photon distribution and spectrum. Using a realistic Galactic stellar distribution model, 
and dust distribution, we calculate the diffuse radiation field from stars in the Galaxy (the 'optical' 
radiation field), including absorption and scattering. Using a dust heating code, we self-consistently 
calculate the infra-red radiation field for the same dust model used for the optical calculation; both 
transient and equilibrium heating are included. We present the calculated radiation field spectra and 
distributions, and will use these to calculate the expected Galactic diffuse IC gamma-ray spectrum. 

1. Introduction 

The Galactic interstellar radiation field (ISRF) plays an important role in the propagation (energy 
losses) and 7-ray emission of cosmic-ray electrons, via the inverse Compton effect. For 7-ray calcu- 
lations we require the spectrum of the ISRF in three dimensions over the whole Galaxy. The most 
detailed calculation to date (Strong et al. 2000) has been widely used, but subsequent to this work a 
large amount of relevant new astronomical information on stellar populations. Galactic structure and 
interstellar dust has become available. Here we give a status report on our efforts to compute a new 
ISRF, exploiting this new information, for use with upcoming 7-ray missions like GLAST. 

2. Modelling the Interstellar Radiation Field 

We use a model for the distribution of stars within the Galaxy based on the statistical SKY model 
of Wainscoat et al. (1992). We have also included the improvements to the model given by Cohen 
(1993, 1994, 1995). The stellar model assumes a type classification, which includes 87 stellar classes 
encompassing main sequence stars, AGB stars, T-Tauri stars, and planetary nebula. For each stellar 
type, there is an associated local density, scale height, fraction of local density in each of several discrete 
spatial components (see later), and spectrum given in the Johnson-Cousins-Glass photometric system, 
plus filters at 12 /im and 25 /im. 

We have made several modifications to the original model to better match recent observations and 
theoretical advances in stellar spectral modelling : 



1. The original model included five discrete spatial components : bulge, arms, disc, halo, and ring. 
Recent analysis of 2MASS data (Lopez-Corredoira et al. 2001) indicates the presence of an 
offset bar, distinct from the bulge. The bar is a thin cylindrical structure with the near side 
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offset with respect to the galactic centre (GC)-Sun Une by I ~ —20° to —30°. We model the bar 
as a thin ellipsoid with parameters as described by Preudenreich (1998), and assume the bar 
stellar population is similar to the disc, with normalisation similar to the bulge. 

2. Analysis of COBE data (Freudenreich 1998) and DENIS star counts (Lopez-Corredoira et al. 
2002) reveals a cut-off within a few kpc from the GC, and no cut-off out to a distance of at 
least 15 kpc; see also Lopez-Corredoira et al. (2004). Ojha (2001) obtains from 2MASS data a 
disc scale-length of ^ 2.8 kpc. Therefore, we model the disc with an exponential cut-off inside 
3 kpc, and a smooth exponential distribution with scale-length 3 kpc for distances greater than 
3 kpc from the GC. 

3. The original SKY model spectra were empirically derived from observations of stars corre- 
sponding to each respective stellar type, or, for AGB stars, from a stellar evolution model. For 
normal stars, we have instead obtained spectra from the synthetic spectral library of Girardi et 
al. (2002). For each stellar type, we obtained reference values for effective temperature, lumi- 
nosity, and surface gravity from Allen (2000). Then, a representative spectrum for each stellar 
type is obtained by interpolating over the Girardi grid. If an appropriate spectrum cannot be 
found in the Girardi grid, the original SKY model spectrum is used. 

The light emitted by stars is absorbed and scattered by dust in the interstellar medium (ISM). 
Absorbed light is re-emitted in the infra-red, while scattered light undergoes further absorption and 
scattering. To calculate the dust extinction (absorption and scattering), and diffuse infra-red emission, 
we require a model for dust in the ISM. We include in our model graphite grains, polycyclic aromatic 
hydrocarbons (PAHs), and silicate grains, since some mixture including these, or similar molecules, 
reproduces the galactic extinction curve (Zubko et al. 2004). We assume that grains in our model are 
spherical, and take the optical absorption and scattering efficiencies for graphite, PAHs, and silicate 
grains from Li & Draine (2001). We adopt the grain model abundance and size distribution from 
Weingartner & Draine (2001) (their best fit Milky Way model). 

Heating of dust grains in the ISM by the ambient optical radiation field produces the diffuse infra-red 
emission. The dust absorption and re-emission process is dependent on the size of the dust grains 
involved. Grains with sizes less than ~ 0.1 /im ('small') are stochastically heated, and undergo 
temperature fluctuations. Grains with sizes larger than 0.1 /im ('large') maintain a nearly constant 
temperature under heating. For 'small' grains, we calculate the transient heating process and re- 
emission using the 'thermal continuous' approach of Draine & Li (2001). For 'large' grains, we obtain 
the infra-red emission by balancing absorption with re-emission as described by Li & Draine (2001). 

Scattering by dust in the ISM produces the so-called diffuse galactic light (DGL) that comprises 
about 10% — 30% of the total optical radiation field (Leinert et al. 1998). We model the scattering by 
assuming a Henyey-Greenstein phase function (Henyey & Greenstein 1941), ^{6,g), for the angular 
distribution of scattered light. The phase function is parameterised by the asymmetry parameter, 
g = g(A), with A the wavelength. We calculate g by summing over grain types, and averaging, 
for each grain type i, the scattering phase function gi{X,a) (a = grain size), which we take from 
Laor & Draine (1993), over our assumed grain size distribution, nj(o), weighted by the scattering 
cross-section, af'''^*^{X,a). 

We assume the dust follows the Galactic gas distribution, and use the gas model for neutral and 
molecular hydrogen described by Strong et al. (2000). Furthermore, the Galactic metallicity gradient 

is assumed to be exponential with radial scale length ^ 4 kpc. 

For the radiation field calculation, we assume a cylindrical geometry. We simplify our calculations 
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by adopting symmetry about the Galactic plane and in azimuth. The maximum radial extent of the 
Galactic volume is taken to be Rmax — 20 kpc, and the maximum height above the plane z^ax is set 
equal to 5 kpc. We divide the total Galactic volume into elements Vi of approximately equal size. For 
each Vi, the radiation field, u(A), with full angular distribution, is calculated as we describe below. 

The optical radiation field is obtained following the method of Kylafis & Bahcall (1987). We first 
calculate the optical radiation field for absorption only (no scattering). This radiation field is then 
used as input to calculate the amount of light scattered only once. The 'once-scattered' radiation 
field is used as input to calculate the amount of light scattered twice, and so forth for the desired 
number of scatterings in the calculation. We have found absorption with the contribution by once- 
and twice-scattered light is sufficient to adequately calculate the optical radiation field. The total 
optical radiation field is obtained by summing these contributions. 

The infra-red radiation field is obtained by using the total optical radiation field for each volume 
element to calculate the emissivity for transient and equilibrium heating. Then, for each Vi, we 
integrate over all volume elements to obtain the infra-red radiation field. Subsequently, the infra-red 
radiation field for each volume element is used to calculate the re-absorbed infra-red emissivity, as 
the optical was used for the infra-red emissivity calculation. Another volume integration is again 
performed for each Vi to obtain the re-absorbed infra-red radiation field. The total infra-red radiation 
field is the sum of these two components. 

Figurenshows our calculated local radiation field, including the contribution by the cosmic microwave 
background (CMBR). Also shown in the figure are observations by Apollo (Henry, Anderson, & Fastie 
1980), DIRBE (Arendt et al. 1998), and FIRAS (Finkbeiner et al. 1999). The agreement with the 
observations is generally good. The scattered component of the optical radiation field is ~ 10 — 20% 
of the total optical radiation field between ~ 0.2 — 2 /im - in good agreement with observations of 
the DGL (Leinert et al. 1998). We note that our calculation falls significantly below the measured 
radiation field in the ultraviolet regime, which may indicate a requirement to include further source 
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Figure 1. Local radiation field modelled as de- 
scribed in the text. Black line : total radiation field, 
including CMBR. Blue solid line : total optical. Blue 
dashed line : total scattered light. Red line : total 
infra-red. Data : thick dot-dashed line, Apollo; thick 
solid line, DIRBE; thick dashed line, FIRAS. 
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Figure 2. Spatial variation of the total radiation 
field throughout the Galaxy. Black line : total radia- 
tion field for GC. Magenta line : total radiation field 
for J? = kpc, z = 5 kpc. Blue line : total radiation 
field for R = 4 kpc, z — kpc. Red line : total ra- 
diation field for R = 12 kpc, z = kpc. Green line : 
total radiation field for R — 20 kpc, z — kpc. 
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classes in the stellar distribution model. Furthermore, the model underpredicts the infra-red from 
~ 10 /im to ~ 40 /im. However, the contamination by zodiacal light around ~ 10 /um in the DIRBE 
observations is uncertain, so this may not be such a problem. 

Figure 13 shows the ISRF at different locations in the Galaxy. The radiation field is most intense in the 
inner Galaxy, and remains so even for distances above the Galactic plane of ~ 5 kpc, and probably 
more. Interestingly, the emission longward of 100 /xm for the in-plane R — Q kpc and R = A kpc 
curves is about the same. This probably reflects the peak in the gas/dust distribution around i? ~ 4 
kpc. However, further analysis will be required to ascertain if this is indeed the true cause of this 
feature. 

3. Diffuse Gamma-Rays 

We will use our new determination for the Galactic ISRF to recalculate the contribution of IC scattered 
gamma-rays to the total diffuse Galactic 7-rays, using the GALPROP cosmic-ray propagation code 
(Strong et al. 2000, 2004). Results will be presented at this conference. 

When completed the new ISRF will be made publicly available as part of the GALPROP package 
(see ^ht t p : / / www . mp e . mpg .de/^aws/aws.htmlj . 
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